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Abstract 

Using the formal analysis made by Bohm in his book, "Quantum theory", Dover Publications Inc. New 
York (1979), to calculate approximately the phase time for a transmitted and the reflected wave packets 
through a potential barrier, we calculate the phase time for a semiconductor system formed by different 
mesoscopic layers. The transmitted and the reflected wave packets are analyzed and the applicability of 
this procedure, based on the stationary phase of a wave packet, is considered in different conditions. For 
the applicability of the stationary phase method an expression is obtained in the case of the transmitted 
wave depending only on the derivatives of the phase, up to third order. This condition indicates whether the 
parameters of the system allow to define the wave packet by its leading term. The case of a multiple barrier 
systems is shown as an illustration of the results. This formalism includes the use of the Transfer Matrix to 
describe the central stratum, whether it is formed by one layer (the single barrier case), or two barriers and 
an inner well (the DBRT system), but one can assume that this stratum can be comprise of any number or 
any kind of semiconductor layers. 
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I. INTRODUCTION 



In the last century the calculation of the time spent by a particle when passing through a potential 
barrier was, for a long time, one of the basic and controversial problems since the early days of 
Quantum Mechanics. When the issue of the delay time of a transmitted wave packet through a po- 
tential barrier was under investigation by MacCollyJ] and later by Hartman,[2] using the Wigner's 
phase time introduced in nuclear physics, the striking superluminal effect arose immediately. At 
this time the question "how much time does tunnellingtake" was loosely formulated. [1] Early 
answers to this problem J3, 4] and alternative proposals 13,^,5,3,7,8,9] run from pure semiclas- 
sical to fully quantum mechanical models. Nowadays the impressive number of low-dimensional 
semiconductors devices brought a new urgency to the essential measurement and/or modelling of 
tunnelling time for charge carriers motion. The last can be seen reflected in the large presence of 
publications. 

Early in the '90s, real experiments on photon-twins interference and on optical pulses 



propagation 



ion[Q, 



111] measured the delay time, in a simple and direct way, at first. On the other 



hand, most of the available experimental setups, pretending to be relevant to the tunnelling is- 
sue, actually involve other times derived from scape and/or decay phenomena. In this sense, 
their results are not able to identify real tunnelling time scale, and consequently should be 
questionable as potentially misleading. Authentically connected to the tunnelling process delay 
measurements. Jlo L [ill fl2fl and uncommonly good agreement with some of them. Jlpj IT] found 
within the phase-time model. fl3| are striking developments from the days of the lively debate on 



these matters appeared during the late 1980's. Approximate [T4Q and multiba nd 11511 pha se time 



1211 and in 



calculations in different systems, confirm experimental results reported in Ref. 
Ref.|Q], respectively. The robustness of the phase time approximation was assured by its con- 
sistency with the Maxwell's equations predictions . 11 1711 Being largely stimulated by the success of 
the phase-time conception, we had applied here the Stationary Phase Method (SPM),[18] firstly to 
evaluate the phase time as it straightforwardly deals with both initial (incident) and final (transmit- 
ted and reflected) dispersion phase amplitudes and finally to further study this magnitude to know 
its applicability more closely. 



[3H 



21, 



22] 



In the last years some other authors have studied these problems in many respects. [1 
At the same time, in Ref. [19] several references of different applications of the tunnelling process 
in different low dimensional semiconductor devices are given. This is the reason to further study 
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this magnitude to know its applicability more closely. 

In the present paper the formal analysis of Bohmyl] is used to determine the transmitted wave 
packet and, properly, the phase time in an arbitrary semiconductor layered system for a charge 
carrier trespassing the structure. The point is that Bohm got the leading term of the transmitted 
wave packet and from it, got an expression to calculate the phase time, and these approximations 
not always are good. To reach when they are good and when do not is our task. The system could 
be described by a single band or multiband models. In figure Q] the general view of the system 
under study is depicted, no matter how many layers are included in the group M. 



Using the Transfer Matrix (TM) method, |23|] the wave packet reflected and transmitted are ob- 
tained, using the SPM II18n to solving the integrals for these waves. This method lead us to an 
applicability condition for it and, properly, the phase time as a function of the parameters of 
the system. The application of the resulting expressions for the Schrodinger single band case is 
given and some results obtained for the double barrier resonant tunnelling semiconductor structure 
(DBRT) are given for illustration. Some comments were included to extend these results to the 
case of a system described by N second order differential system. 



H. THE FORMAL ANALYSIS 



In the system depicted in figureCQthe Schrodinger wavefunction can be written as: 



1po(z,PR,t) 



D( PR )" 1 
+F{p R 



^o(z,p R ) 
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where the TM of wavefunction and derivative |24fl was included to describe the central layer 



M.[|23j] As this part can be arbitrary, the expression for the TM will depend on the form of the 
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potential of this layer. Also it was written p R = \/2m*E. The wavefunction for the layer L was 
also written in terms of p R for convenience. In doing this, the potential Vlr was defined (see 
caption of Figure [T]). 

The wave packet is obtained, for different values of coordinate z by forming the expression: 



*{z,t) 



dpRf(pR-pRo)lpo(z,PR,t), 



0) 



where function f(pR — pr ) is a shape function which peaks at the value pr q and rapidly goes to 
zero for large values of the difference p R — p Ro , then the integral limits can be extended to ±00. 
In the case of the transmitted wave, taking D(p R ) = 1 as a condition of normalization of the 
wavefunction used to form the wave packet, one obtains for ^ R (z, t) for z in region R: 

/oo 
dp R f{pR - Pro)B( P r) Q{^p R m+*M PR )-i(E( PR ) )m (4) 
-00 

The normalization condition means that the incident wave to form the packet is normalized in the 
L region. Here it was written A(p R ) = B(p R ) Q^pr) using 4>t(p R ) as the phase of the transmitted 
wave amplitude. For the case N > 2 (i. e. a physical system described by two or more coupled 
differential equations), the condition of normalization have to be released because it is necessary 
to write the spinor as a part of the wave function. j25\\ In parameter A(p R ) the matching of the 
different layers in the structure is included. For N > 2 this parameter is a vector, then this 
matching process appears in the coefficient B(p R ) and in the phase 4>t{pR) which in the multiband 
case must be calculated by components and no matrix expression can be given. 
Considering the SPM to perform the integral in the case N — 1 one has to expand in Taylor's 
series the exponent in © (which we called Qt(pr)) and taking the value of p R which produces an 
extreme for the exponent ipno) as the approximation, one obtains: 



®t(pr) 
<£>t(pro) 
V R {z,t) 



G{p 



Ro 



i(zp R /h) + i(/)t(p R ) - i {E(p R ))t/h. 

The definition of the exponent 
i{zp Ro /K) - i {E(p Ro )) t/h + i(p t (PRo) + 0(p 

The Taylor series up to second order, 

G(p Ro ) Q^-^PRo)/ 11 ) '+( z PRo/H)] gi(<MPfl )+7r/4) 

The leading term of the wave packet, 

V^B{p Ro ) 



L) + 



(d 2 (j) t (p R )/dp R ] 



PRo 



[l/(m*v gR )] {d(j) t {p R )/dp 
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(5) 



(6) 



(7) 



(8) 



The coefficient of the leading term of the wavepacket. 



Here we have considered that f(pR — Pro) = 1 and we use v 9 r as the group velocity of the 
packet in layer R. Expression © is the leading term of the transmitted wavefunction, obtained by 
making this approximation. The coefficient of this wavefunction is given by ([8]) written in terms 
of the phase delay of the transmitted wave. 

The applicability of the SPM takes into account that it uses the Taylor expand of the exponent and 
neglects the terms from the second order. This leads to write: 

1 



<3>t(pr) 
®t(Pro)(pr -Pro) 3 



&t(pro) + &'t(Pro)(pr ~ Pro) + ^®t(PRo)(Pr ~ Pro? + ■ ■ • • (9) 



Qt(pr ) + -®t(Pro){Pr - Pro) 2 



®t(Pro)(pr- Pro) 



< e; 



with e — > 0. 



as is an extreme .(10) 
(ID 



<3>'t(PRo) (pr ~ Pro) 2 <S>t(pro) 
Evaluating the derivatives of the exponential phase (fTTI) in terms of the derivatives of the phase of 
the transmitted wave one obtains as the condition for the applicability of the SPM the expression: 



dpi 



1 / dchipR, 

PRo 



dpR 



+ 
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<C 1. (12) 



Pr=PRo \ ' Pr=PRo 

This is the main contribution of this paper. This condition evaluates the applicability of the SPM 
and points over the use of the phase delay time for every group of values of the parameters of the 
system. 

Nevertheless, this expression (fT2j) has the numerator and the denominator dimensional and the 
quotient non dimensional, then to properly compare these expressions it is better to multiply by 
p Ro both, numerator and denominator. This lead us to: 
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(13) 
(14) 

(15) 



The phase time for the transmitted wave is obtained from the condition of stationary phase of the 
exponential in the integral ©. After including the matching at layer boundaries, one has for the 
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phase of the exponential in © the expression: 



ar(PR,t) = + Mpr) ~ t e (Pr), (16) 

ft n 

daT - =► r T = ^ ( ^M] ; (17) 



d PR PRo V dp„ /pR=pRo 

which is the formula to evaluate the phase delay time of the transmitted wave. [3] In (fT6l ) was taken, 
as Bohm did in his book, z = Az = because it refers to the phase between group M of layers 
and layer R, i. e., the wave packet reaches the same position, later than if there were no dispersion 
potential causing the wave to be reflected. In this sense, the phase of layer R differentiates, bearing 
a term that comprises wave packets evolution delay information. For the case of N > 2 bands, 
the whole analysis cannot be generalized for the present scheme from the case N = 1 because 
the step of converting a complex number from a + ib to pQ 1 ^ cannot be performed in matrix 
notation and one must passes to 2N x 2N components. Further investigation is required to write 
close expressions in this case. This is important because there are several problems described by 



the standard Sturm-Liuoville N x N differential equation system ll23ll of gr eat p ractical interest. 
Models as that due to Bogoliubov for superconductor excitations description ||26| l could be treated 
as well. 

A simple consideration of closeness between the phase-time model and the dwell time (within its 



phase-time probabilistic average formulation Ii27l l28h). dispose us to speculate that the requirement 



(1121) should be readily suited to it, with minor changes. We are interested in comparing these two 
possible conditions to get light into the use of different times for tunnelling processes. 

HI. RESULTS AND DISCUSSION 

The application of this formal analysis to different physical systems allows one to determine 
whether the phase time can be applied to a given system and to obtain it from the wavefunc- 
tion. As an illustration we applied this procedure to the case of a double barrier resonant structure 
device in GaAs/AlGaAs considering the parameters shown in Table HI 

The potential of the system is depicted in Figure [H where the extreme left and right layers were 
considered as metallized contacts, which are semiconductors (GaAs) with flat band and an electric 
field applied to the structure. 

Using (flTl) . after making the matching considering the differences of masses in each layer by using 



the TM algorithm. Il23ll the phase delay time has the behavior depicted in figure [3] as a function of 



the energy of the incident wave. 

Our results for the phase time depicted in Figure [3] are of the same order of magnitude of other 
calculations and the behavior of the phase time is as others achieved, as can be seen in table fll] for 



electrons and photons in similar system, reported elsewhere. 
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29 



30 



3 IN Several methods were 



time,|29j] Wentzel- Kramer- Brillouin (WKB) 



used by these authors, namely: lifetime, [29j] dwel 
quasi-classical approximation |30J] and phase time. II 1311 In the case of photons, the reported values 
correspond to 1.5 /im optical pulse wavelength, propagating through double-barrier photonic band 
gap (FBG) J3 in In this table are included some useful data as if there is applied electric field, if the 
results were achieved theoretical or experimentally and the model used to perform the calculation. 
The application of the SPM to this system is governed by expression (PT2l) and in Figures |4] a, b, 
c, d and e are shown separately the numerator Q N (p R )/l0 6 , the denominator Q D (p R ) /10 14 and 
the quotient XQor) of the applicability condition (fT3l) for different energy ranges. This is the main 
result of this paper, because the phase delay time approximation is already known, but (fT3l) is not 
used to assure its application to different systems. It is easily seen that in all graphics is under 
Qd, so the procedure and the phase delay time are valid for all the energy range of interest. Never- 
theless, there is an isolated point, seen in figure 4e) that goes over unity, and makes the SPM and 
the phase delay time inapplicable. 

This analysis allows to say that this definition of time is good enough for many useful analysis at 
all energy ranges. 

As a conclusion we have calculated the phase delay time in a system of semiconductor layers, 
illustrating with the simple case of a DBRT system described by the Schrodinger equation with an 
electric field applied and some light about the applicability of this definition of time is given by 
considering the condition obtained for the use of the SPM in reaching the transmitted wave packet. 
It is clear that one has to apply the applicability condition in each case under study to assure that 
the phase delay time is good in the conditions of each concrete problem. It is also an interesting 
guess to extrapolate the applicability condition obtain ed (fl2l) for the phase time, to the case of the 
dwell time in its probabilistic average formulation ||27[ 128( 1 with minor changes. 
Also in this paper some considerations were made to extend these formulae to the case of systems 
with N second order coupled differential equations in which some of the algebra must be done in 
matrix notation and other cannot. This application leads to individual results for each component 
separately and after that one can rebuild the matrices. Expression (fT2~l) is valid for each one of the 
components and must be obtained and evaluated individually. The application of these results to 
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the case of N second order differential equations is in progress. 
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TABLE I: Parameters of the DBRT considered in the calculation sketched in Figure [2] 



No Parameter 


Value 


1 


Barrier Height (Vi) 


250 meV 


2 


Difference between the band ed| 


^es 




of sides L and R (V Q ) 


40meV 


3 


Barrier width b = z\ — a\ 


40 A 


4 


Well Width d = z 2 - z x 


100 A 


5 


m\ in units of m 


0.066 


6 


m\ in units of m 


0.8 
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TABLE II: Time scale for tunneling 



System 


Potential Structure 


Data Source 


Resonance 


Bias 


Time 


Value of Time 










[eV] 


[meV] 




[ps] 


electrons 


DBRT 


Alo.3Gao.7As/GaAs 


Teo. Ref[29] 


0.05 


- 


life 


5.7 


electrons 


DBRT 


Alo.3Gao.7As/GaAs 


Teo. Ref[29] 


0.18 


- 


life 


2.9 


electrons 


DBRT 


Alo.3Gao.7As/GaAs 


Teo. Ref[29J 


0.05 


- 


dwell 


5.7 


electrons 


DBRT 


Alo.3Gao.7As/As 


Teo. Ref[29] 


0.18 




dwell 


2.9 


electrons 


DBRT 


Gao.47Ino.53As/Alo.48lno.52 


As Teo. Ref[30] 




40 


WKB 


0.5 


electrons 


DBRT 


Alo.3Gao.7As/GaAs 


Teo. Ref[13] 


0.11 




phase 


0.02 


electrons 


DBRT 


Alo.3Gao.7As/GaAs 


Teo. Ref[13] 






phase 


ps 0.02 


photons 


FBG 


mono-mode optical fiber 


Exp. Ref[3TJ 






traversal 


180 


photons 


FBG 


mono-mode optical fiber 


Teo. RefQJJ 






phase 


300 
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Figure Captions 
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FIG. 1: General view of the Potential System under study. Layers L and i? must not be equal necessarily. 
Vlr is the potential difference between the L band edge and the and the R band edge taken as energy 
reference level. In z-axis the interfaces are named as ax and 0,2- 




0, Z< Z 2 Q 2 

FIG. 2: General view of the DBRT system under study. Here V Q is called the Vlr of FigureQ] Points a\ and 
a 2 are here the same as in Figure [TJ then points z\ and 22 and potential V\ belong to layer M of FigurdT] 
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Vl=250 meV 
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FIG. 3: Phase delay time calculated by (fTTT ) for the DBRT with the parameters given in Table U For E < V Q 
the picks correspond to energies of the inner well. 
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FIG. 4: Evaluation of the condition of applicability of the SPM (PT2l)in the energy range. The shown graphs 
are: a) Denominator (Qd{pr) /10 14 ) versus energy, b) The same as a) but in a particular energy range, c) 
Numerator (Qn(pr)/10 6 ) vs energy, d) The same as c) but in a smaller energy range, and e) Quotient 
versus energy in the whole energy range. 
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